Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1278 1611 (N41W25 L=277)   305     45        0 (1256/1587-1299/1636)  38.7 W22-W27 N39-N43
      2 1784 1299 (N16W60 L=312)  1048    140        2 (1722/1254-1835/1345) -35.0 W54-W65 N12-N19
      3 1843 1222 (N11W65 L=317)   487     45        1 (1827/1190-1866/1247)  53.1 W61-W68 N08-N12
      4  223  936 (N00E60 L=192)   384     55        0 ( 205/ 905- 249/ 963)  26.6 E62-E56 S01-N02
      5 1154  867 (S06W07 L=259)   253     97        0 (1108/ 837-1193/ 888) -31.0 W04-W09 S07-S04
      6  221  659 (S18E69 L=183)   789     56        1 ( 198/ 635- 251/ 684) -45.0 E73-E63 S19-S16
      7 1416  583 (S26W25 L=277)  1280    309        2 (1296/ 488-1518/ 678)  40.6 W17-W32 S32-S20
      8  627  588 (S22E29 L=222)   412    125        0 ( 592/ 550- 682/ 626) -45.0 E31-E25 S24-S19
      9 1041  519 (S28E01 L=250)   311    137        0 ( 989/ 506-1107/ 555)  23.2 E04-W03 S29-S25
     10  590  517 (S27E34 L=218)   297    104        0 ( 566/ 474- 628/ 549)  45.0 E36-E30 S30-S24
     11  770  487 (S30E21 L=231)   352    119        1 ( 754/ 433- 785/ 542)  76.0 E22-E19 S33-S25
     12  929  468 (S32E10 L=242)   318     66        0 ( 901/ 442- 955/ 498)  53.1 E12-E07 S33-S29
  

    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