Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  919 1787 (N60W14 L=287)   381     48        0 ( 900/1759- 937/1812)  51.3 W11-W18 N57-N62
      2 1530 1612 (N31W53 L=326)  2197    265        2 (1425/1579-1660/1644) -11.8 W42-W66 N25-N34
      3  605 1469 (N42E23 L=249)  1059    226        0 ( 506/1449- 716/1489) -11.9 E32-E13 N39-N44
      4  248 1456 (N45E69 L=204)  1224     97        1 ( 204/1423- 302/1474)   0.0 E85-E58 N42-N46
      5 1767  676 (S31W49 L=322)  2987    378        4 (1653/ 560-1878/ 843)  35.0 W37-W62 S37-S22
      6  234  629 (S07E68 L=204)  1267     87        3 ( 202/ 584- 260/ 660) -51.3 E72-E63 S09-S04
      7  398  526 (S15E54 L=218)   488    120        0 ( 366/ 486- 444/ 567) -56.3 E56-E51 S17-S11
      8  656  469 (S21E36 L=237)   273     40        0 ( 635/ 444- 680/ 493) -76.0 E36-E34 S23-S19
      9 1077  279 (S44E14 L=259)  3232    489        3 ( 871/ 189-1296/ 394)   9.5 E36-W06 S48-S38
  

    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