Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  432 1491 (N31E49 L=179)   334     49        1 ( 411/1469- 460/1517) -15.9 E53-E46 N29-N32
      2  541 1193 (N12E32 L=196)   233     65        0 ( 526/1159- 557/1230) -63.4 E33-E31 N10-N14
      3 1775  953 (S01W54 L=282)   405     78        0 (1743/ 927-1802/ 985) -26.6 W50-W57 S02-N01
      4  428  705 (S19E42 L=187)   364    117        0 ( 377/ 690- 482/ 721)   0.0 E46-E36 S20-S17
      5  199  667 (S23E72 L=156)   853     66        2 ( 176/ 642- 226/ 698)  15.9 E82-E66 S25-S21
      6 1706  655 (S20W52 L=281)   351     50        0 (1685/ 635-1728/ 678)  56.3 W50-W54 S21-S18
      7  435  611 (S25E43 L=185)   525    130        0 ( 417/ 552- 470/ 664)  53.1 E46-E38 S29-S21
      8  811  511 (S31E14 L=214)   248     60        0 ( 776/ 499- 844/ 524)   0.0 E16-E11 S31-S29
      9 1522  457 (S34W42 L=270)   322     15        0 (1501/ 437-1541/ 477)  18.4 W40-W44 S35-S32
     10 1147  368 (S41W12 L=240)   515     56        2 (1117/ 343-1179/ 392)  31.0 W09-W14 S43-S39
     11  896  326 (S46E09 L=219)   900    168        2 ( 815/ 286- 954/ 358) -22.6 E16-E04 S49-S43
  

    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