Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1027 1526 (N28W00 L=  5)   388     75        0 ( 983/1489-1077/1551) -38.7 E02-W03 N25-N29
      2  268 1345 (N18E54 L=311)   696    108        0 ( 253/1298- 283/1401) -66.8 E56-E52 N14-N21
      3 1686 1343 (N17W45 L= 51)   489     91        1 (1645/1320-1723/1360) -18.4 W42-W48 N15-N18
      4  163 1211 (N10E65 L=301)   742    161        1 ( 122/1192- 202/1226)   0.0 E71-E59 N08-N10
      5 1790 1138 (N05W53 L= 58)   662    132        2 (1749/1102-1832/1173) -20.6 W49-W57 N02-N06
      6  650  885 (S11E23 L=342)   689    203        0 ( 617/ 781- 684/ 987)  67.4 E26-E20 S17-S05
      7 1071  973 (S06W03 L=  8)   183     39        0 (1044/ 959-1099/ 987)   0.0 W01-W04 S07-S05
      8  207  697 (S21E63 L=303)  1994    427        2 ( 187/ 582- 243/ 845)  53.7 E68-E57 S27-S12
      9 1010  685 (S24E01 L=  4)   403     63        1 ( 987/ 642-1037/ 735) -80.5 E02-W00 S26-S20
     10 1846  683 (S22W63 L= 68)   573    112        3 (1809/ 652-1882/ 723)  32.0 W59-W67 S24-S19
     11  417  631 (S27E44 L=321)   919    204        2 ( 389/ 558- 449/ 718)  55.0 E47-E40 S31-S21
  

    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