Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  622 1669 (N22E44 L=238)  1684    187        5 ( 541/1619- 718/1703)  35.0 E49-E38 N17-N25
      2 1533 1507 (N35W20 L=303)  1268    227        1 (1419/1448-1623/1558)   0.0 W11-W28 N33-N37
      3  288 1501 (N06E65 L=218)  1018    103        6 ( 243/1457- 320/1545)  51.3 E69-E61 N02-N08
      4 1797 1208 (N25W47 L=329)  1131    259        1 (1767/1104-1822/1340) -86.2 W43-W50 N17-N32
      5  192 1270 (S09E65 L=217)   646     97        1 ( 174/1230- 220/1317) -32.0 E69-E60 S11-S06
      6  813 1221 (S01E16 L=266)   357     92        2 ( 763/1192- 851/1244)  11.3 E19-E14 S02-N00
      7 1719 1144 (N18W40 L=322)   277     33        1 (1698/1127-1745/1165)  14.0 W38-W42 N16-N19
      8 1847  815 (N05W61 L=343)   495    112        0 (1799/ 792-1888/ 854)  36.9 W55-W64 N02-N08
      9  987  755 (S23W05 L=287)   281    100        0 ( 960/ 726-1030/ 789) -21.8 W02-W07 S24-S21
     10  489  574 (S47E28 L=254)   798    137        1 ( 431/ 537- 555/ 610)  -4.4 E34-E21 S48-S45
     11  625  540 (S46E15 L=267)   361     67        2 ( 597/ 512- 651/ 573) -45.0 E17-E12 S47-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