Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  787 1530 (N33E11 L=108)   446    101        1 ( 741/1488- 824/1565) -39.8 E14-E08 N29-N35
      2  160 1298 (N25E63 L= 56)  1248    195        1 ( 134/1236- 176/1376)  90.0 E65-E60 N21-N30
      3 1121 1213 (N09W08 L=126)   251     89        0 (1077/1178-1162/1244) -51.3 W05-W09 N06-N11
      4  481 1154 (N12E33 L= 86)   288     12        0 ( 462/1129- 500/1177)  71.6 E34-E31 N10-N13
      5 1887  850 (S19W62 L=181)   463    103        1 (1853/ 831-1912/ 874)   8.1 W58-W65 S20-S17
      6 1624  742 (S24W38 L=156)   331    160        0 (1568/ 721-1687/ 769) -12.5 W33-W42 S24-S21
      7  296  710 (S12E54 L= 64)   651    101        2 ( 263/ 673- 335/ 744) -23.2 E57-E50 S14-S10
      8 1092  511 (S33E01 L=118)   327    114        0 (1070/ 459-1108/ 569)  66.8 E03-W00 S37-S29
      9 1464  487 (S39W27 L=146)   852    204        3 (1379/ 437-1549/ 555)  26.6 W20-W33 S42-S35
     10  307  473 (S27E63 L= 56)  1375    253        2 ( 273/ 406- 333/ 537)  56.3 E67-E59 S31-S22
  

    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