Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  720 1297 (N23E19 L=195)   299    100        2 ( 696/1266- 744/1330) -53.1 E20-E17 N21-N25
      2  768 1180 (N16E16 L=199)   250    111        0 ( 715/1160- 821/1199)  -8.1 E19-E12 N14-N16
      3  706  933 (N01E20 L=194)   280     42        3 ( 683/ 911- 733/ 957) -18.4 E21-E18 N00-N02
      4 1226  891 (S05W12 L=226)   620    188        2 (1142/ 859-1307/ 932) -26.6 W06-W16 S07-S01
      5  198  778 (S09E67 L=148)  1236    122        2 ( 165/ 735- 232/ 829) -50.2 E71-E62 S11-S05
      6 1636  676 (S22W42 L=256)   445    127        0 (1598/ 622-1664/ 736)  66.8 W39-W44 S25-S18
      7  636  608 (S20E29 L=185)   638    158        0 ( 595/ 554- 700/ 655)  26.6 E33-E23 S23-S16
      8  371  604 (S20E52 L=163)   490     83        1 ( 348/ 567- 397/ 639) -90.0 E53-E49 S22-S17
      9  553  442 (S32E41 L=173)  1382    247        2 ( 530/ 316- 590/ 543)  54.2 E51-E34 S42-S24
     10 1624  495 (S35W46 L=261)   364     63        0 (1595/ 478-1657/ 510)  18.4 W43-W49 S36-S34
  

    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