Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1648 1452 (N20W52 L=187)   536    138        1 (1598/1427-1698/1469) -18.4 W47-W56 N17-N21
      2 1171 1196 (N14W12 L=146)   424    128        0 (1147/1123-1196/1257)  77.5 W10-W12 N09-N18
      3 1611 1028 (S04W38 L=172)   545    122        1 (1553/1000-1661/1051) -14.0 W33-W41 S05-S01
      4 1646  773 (S20W38 L=172)  1352    347        3 (1554/ 608-1700/ 932)  57.7 W29-W42 S29-S10
      5 1887  862 (S21W62 L=197)   619     97        0 (1867/ 812-1905/ 905)  51.3 W59-W65 S23-S17
      6  281  764 (S01E57 L= 77)   351     75        0 ( 267/ 727- 296/ 798)  53.1 E58-E54 S03-N01
      7 1293  676 (S19W11 L=146)   309    103        1 (1280/ 625-1309/ 731)  74.1 W10-W12 S22-S15
      8  302  690 (S06E57 L= 77)   367     57        1 ( 281/ 671- 333/ 709) -21.8 E59-E53 S06-S04
      9  347  608 (S11E55 L= 79)   433     58        1 ( 328/ 583- 366/ 635) -90.0 E57-E52 S12-S09
     10  784  506 (S22E24 L=110)  1079    243        2 ( 684/ 452- 887/ 559)  13.2 E33-E16 S25-S19
     11  568  321 (S34E50 L= 85)   817    175        1 ( 538/ 290- 592/ 358)  33.7 E53-E44 S36-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