Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  642 1657 (N25E44 L=130)  2980    655        6 ( 403/1541- 864/1719)  28.1 E58-E27 N13-N31
      2  616 1460 (N12E36 L=137)   267     19        0 ( 600/1437- 638/1479)  90.0 E37-E35 N10-N13
      3  231 1296 (S05E63 L=110)   483     91        1 ( 208/1257- 249/1327)  63.4 E65-E61 S07-S03
      4  689 1240 (N01E25 L=149)   351     86        1 ( 649/1213- 730/1269)  45.0 E26-E22 S01-N03
      5 1214 1098 (N05W09 L=182)   278     84        2 (1194/1051-1237/1143) -68.2 W07-W09 N02-N07
      6 1744  714 (S02W57 L=230)  1091    204        4 (1687/ 646-1775/ 789) -22.6 W49-W61 S05-N01
      7 1217  740 (S15W18 L=192)   319    105        1 (1188/ 705-1257/ 767)  36.9 W16-W20 S16-S12
      8 1464  656 (S13W37 L=210)   612    142        0 (1409/ 604-1523/ 719)  71.6 W35-W38 S17-S08
      9 1194  672 (S19W19 L=193)   214     91        0 (1153/ 659-1245/ 683)  31.0 W16-W22 S20-S17
     10  874  497 (S39W05 L=179)  1088    376        1 ( 710/ 472-1074/ 524)  19.8 E07-W18 S42-S32
     11 1587  467 (S20W55 L=228)   494    107        0 (1534/ 431-1632/ 501)  50.2 W52-W57 S23-S16
     12 1326  363 (S34W42 L=216)   383    133        1 (1290/ 341-1370/ 401)  -6.3 W37-W46 S35-S32
  

    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