Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  406 1654 (N24E68 L=231)   733     56        1 ( 380/1633- 432/1679)  -8.1 E77-E63 N22-N25
      2  607 1647 (N27E46 L=253)   306     95        0 ( 566/1622- 639/1663)   0.0 E50-E42 N25-N28
      3  859 1500 (N23E22 L=278)   277    109        0 ( 838/1452- 879/1547) -39.8 E24-E18 N20-N25
      4 1213 1442 (N27W03 L=302)   885    317        1 (1088/1346-1328/1497) -13.2 E06-W11 N23-N29
      5 1390 1421 (N30W15 L=315)   259    106        1 (1357/1379-1426/1466) -33.7 W12-W18 N27-N32
      6 1883 1183 (N27W65 L=  5)  2679    235        2 (1848/1075-1906/1289)  81.9 W61-W71 N20-N34
      7  451 1061 (S11E37 L=262)   526     92        1 ( 402/1033- 490/1086)  33.7 E40-E34 S13-S09
      8 1771  984 (N12W52 L=351)   666     99        0 (1745/ 934-1792/1031)  74.1 W49-W53 N08-N15
      9 1453  547 (S21W40 L=339)  2178    410        4 (1374/ 429-1569/ 645)  82.9 W35-W50 S28-S12
     10  476  551 (S44E32 L=267)   342     35        1 ( 460/ 524- 497/ 579)  68.2 E34-E29 S46-S41
     11  920  527 (S34W05 L=305)   583    196        1 ( 833/ 497-1009/ 555)  28.6 E01-W10 S37-S30
  

    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