Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1140 1627 (N36W24 L=177)  1292    144        4 (1048/1571-1257/1662) -15.9 W16-W32 N32-N37
      2 1781 1381 (N05W61 L=214)   881    111        2 (1741/1336-1822/1427) -60.3 W57-W63 N01-N08
      3  160 1308 (N38E60 L= 92)   691    251        0 ( 122/1281- 198/1324) -18.4 E64-E55 N35-N38
      4 1762 1179 (S05W51 L=204)   365     66        2 (1736/1150-1783/1209)  31.0 W48-W53 S06-S03
      5  457 1142 (N21E33 L=120)   235     91        0 ( 431/1122- 489/1158)   0.0 E35-E30 N20-N22
      6 1322  974 (S06W16 L=169)   928    316        3 (1242/ 807-1386/1130)  45.0 W07-W22 S14-N01
      7 1790 1092 (S11W51 L=204)   383     73        0 (1754/1068-1815/1114) -31.0 W48-W54 S12-S08
      8  494  944 (N09E34 L=119)   229     62        2 ( 473/ 908- 512/ 979)  56.3 E35-E32 N06-N11
      9  931  455 (S29E20 L=133)   692    187        1 ( 903/ 374- 953/ 532)  65.6 E24-E17 S34-S23
     10 1348  476 (S38W08 L=161)   506    102        0 (1299/ 430-1389/ 512)  26.6 W03-W11 S40-S36
     11  550  347 (S28E51 L=101)   909    227        2 ( 528/ 289- 573/ 432)  61.4 E55-E47 S32-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