Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  664 1602 (N44E29 L=261)   308     61        1 ( 638/1576- 698/1623) -26.6 E32-E24 N41-N45
      2 1121 1429 (N30W09 L=299)   589    124        3 (1087/1383-1149/1492)  45.0 W05-W11 N26-N33
      3  290 1415 (N30E62 L=229)   504    131        0 ( 255/1367- 338/1451)  45.0 E65-E57 N26-N32
      4  298 1337 (N25E57 L=233)   447     68        3 ( 280/1305- 318/1371)  78.7 E58-E54 N22-N27
      5 1203 1331 (N23W13 L=304)   484    113        0 (1147/1312-1256/1354) -23.2 W09-W16 N21-N24
      6 1751  854 (S11W51 L=342)   602     43        1 (1724/ 825-1775/ 883) -63.4 W48-W53 S12-S08
      7 1079  810 (S09W03 L=293)   377     99        1 (1050/ 771-1105/ 855) -59.0 W00-W04 S11-S06
      8  894  774 (S11E09 L=281)   297     71        1 ( 862/ 745- 928/ 799) -26.6 E11-E07 S12-S09
      9 1348  747 (S14W20 L=310)   247     70        0 (1334/ 709-1358/ 787)  90.0 W18-W20 S16-S11
     10 1703  624 (S26W51 L=341)  1121    195        1 (1674/ 580-1745/ 692) -58.0 W48-W56 S28-S20
     11  517  652 (S18E37 L=254)   386    103        2 ( 469/ 630- 563/ 676)   7.1 E41-E33 S19-S16
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y