Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  842 1618 (N23E27 L=242)   389    108        0 ( 780/1587- 890/1639)  45.0 E30-E24 N19-N25
      2  356 1581 (N11E62 L=206)   654     64        1 ( 332/1558- 391/1607)  31.0 E66-E58 N09-N12
      3 1657 1444 (N35W32 L=300)  1712    393        1 (1500/1342-1777/1581) -13.0 W16-W44 N32-N40
      4  363 1398 (N00E52 L=217)   994    258        0 ( 288/1317- 445/1506)  63.4 E56-E48 S05-N07
      5  328 1234 (S10E50 L=219)   535     90        1 ( 291/1208- 366/1266)  -7.1 E53-E45 S11-S08
      6  934 1189 (S00E09 L=260)   203     86        0 ( 919/1149- 949/1228) -45.0 E10-E07 S02-N01
      7 1853 1061 (N18W54 L=323)   370     22        0 (1832/1040-1875/1080)  45.0 W52-W57 N16-N19
      8 1188  666 (S23W19 L=288)   608    145        1 (1147/ 633-1249/ 724) -14.0 W15-W23 S24-S20
      9 1160  609 (S27W20 L=288)   229     50        0 (1137/ 582-1180/ 643) -26.6 W17-W21 S27-S25
     10  705  498 (S46E07 L=262)  1560    399        1 ( 555/ 434- 891/ 569)   3.8 E21-W09 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