Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1269 1586 (N39W05 L=148)   270     68        0 (1245/1566-1303/1605)  26.6 W03-W07 N37-N40
      2  938 1469 (N24E15 L=129)   702    289        1 ( 852/1353-1033/1564) -29.4 E22-E06 N18-N28
      3 1385 1529 (N38W15 L=158)   309     69        0 (1365/1504-1417/1546)  45.0 W13-W17 N35-N39
      4 1677 1379 (N34W40 L=184)  2115    368        1 (1546/1310-1799/1458)  -7.4 W28-W52 N30-N37
      5 1399 1174 (N15W20 L=164)  1130    369        2 (1351/1047-1497/1257)  90.0 W16-W26 N07-N20
      6 1147  982 (S01W08 L=152)   458    140        0 (1061/ 953-1204/1022)  35.5 W03-W10 S04-N01
      7  500  708 (S31E30 L=114)   326     75        0 ( 478/ 674- 524/ 743)  59.0 E31-E27 S33-S28
      8 1153  695 (S18W14 L=158)   233     36        1 (1130/ 670-1173/ 716) -26.6 W12-W16 S19-S17
      9 1616  661 (S11W47 L=191)   825    179        3 (1558/ 634-1676/ 694)   0.0 W41-W52 S12-S09
     10 1477  385 (S31W49 L=193)   334     34        0 (1462/ 364-1494/ 408) -45.0 W46-W51 S32-S29
     11 1245  251 (S47W41 L=185)   375     59        1 (1220/ 232-1267/ 275)  14.0 W39-W44 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