Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  922 1470 (N24E16 L=129)   634    249        0 ( 836/1377-1006/1553) -28.1 E23-E08 N19-N27
      2 1369 1534 (N38W13 L=159)   255     43        2 (1348/1512-1395/1553)  26.6 W11-W15 N35-N39
      3 1489 1489 (N37W24 L=169)   293     69        0 (1460/1471-1524/1504)  31.0 W21-W26 N35-N38
      4 1669 1381 (N34W39 L=184)  1952    360        2 (1542/1309-1788/1471) -10.3 W28-W51 N30-N37
      5 1391 1161 (N14W20 L=165)  1166    401        5 (1329/1049-1472/1264) -33.7 W14-W26 N07-N20
      6 1132  995 (S01W07 L=152)   307    111        0 (1067/ 965-1167/1035)  45.0 W03-W08 S03-N01
      7  470  695 (S32E32 L=113)   516    165        0 ( 426/ 639- 507/ 756)  56.3 E36-E29 S36-S27
      8 1131  703 (S18W13 L=158)   263     34        0 (1111/ 682-1156/ 728) -33.7 W11-W14 S19-S16
      9 1596  668 (S11W45 L=190)   761    187        1 (1532/ 638-1662/ 700)   4.8 W39-W51 S12-S09
     10 1450  397 (S31W46 L=191)   686    177        2 (1392/ 370-1502/ 430)   0.0 W40-W51 S32-S29
     11 1234  258 (S46W40 L=185)   367     57        1 (1209/ 241-1258/ 282)  76.0 W37-W42 S48-S44
  

    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