Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1009 1710 (N51E03 L=212)   351     50        0 ( 983/1687-1040/1730)  31.0 E05-W00 N48-N52
      2  565 1460 (N30E36 L=179)   696    151        3 ( 509/1401- 607/1525) -32.5 E42-E30 N25-N34
      3  154 1264 (N14E77 L=138)   903     89        2 ( 135/1217- 167/1308) -71.6 E81-E73 N11-N17
      4  611  692 (S19E27 L=187)   267     99        0 ( 562/ 677- 656/ 709)   7.1 E31-E23 S19-S17
      5  289  629 (S25E59 L=156)  1559    164        5 ( 253/ 549- 326/ 702)  35.2 E71-E52 S31-S19
      6  605  605 (S25E29 L=186)   440    162        2 ( 569/ 535- 660/ 651)  53.1 E32-E24 S29-S21
      7 1574  600 (S24W41 L=256)   328    101        0 (1540/ 558-1601/ 644)  63.4 W39-W42 S26-S20
      8 1302  358 (S42W25 L=240)   549     76        0 (1275/ 335-1326/ 395)  78.7 W22-W28 S44-S38
      9 1048  318 (S46W04 L=218)  1059    230        0 ( 952/ 283-1146/ 360) -17.4 E04-W12 S49-S42
  

    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