Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1121 1723 (N52W08 L=211)   351     59        1 (1097/1695-1160/1742)  18.4 W06-W12 N49-N53
      2  548 1517 (N34E39 L=163)  1810    406        2 ( 398/1445- 723/1599) -11.6 E61-E22 N29-N39
      3  224 1247 (N14E63 L=139)   664    107        1 ( 210/1195- 245/1293) -80.5 E65-E61 N11-N17
      4  734  687 (S18E19 L=184)   504    215        0 ( 637/ 664- 827/ 721)  -8.7 E25-E12 S19-S16
      5  397  586 (S27E48 L=154)  2214    308        3 ( 319/ 439- 469/ 690)  32.0 E73-E38 S39-S19
      6 1703  645 (S21W52 L=254)   351     42        0 (1687/ 626-1722/ 668) -45.0 W50-W54 S22-S19
      7  794  610 (S23E15 L=187)   448    175        0 ( 747/ 550- 871/ 651)  37.9 E18-E09 S27-S20
      8 1176  497 (S31W12 L=214)   267     83        0 (1135/ 479-1219/ 515)  18.4 W08-W14 S32-S29
      9 1569  413 (S38W49 L=252)   791    114        0 (1538/ 370-1591/ 469) -63.4 W46-W52 S41-S33
     10 1456  363 (S42W40 L=242)  1106     93        3 (1408/ 336-1517/ 396)   0.0 W34-W46 S44-S39
     11 1216  308 (S47W19 L=221)  1207    199        3 (1124/ 285-1310/ 347)  -3.4 W10-W27 S48-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