Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1775 1505 (N16W66 L=188)  1755    310        1 (1695/1435-1858/1544) -29.4 W57-W74 N09-N18
      2 1202 1285 (N13W15 L=137)   427    166        1 (1133/1238-1296/1326) -41.2 W11-W19 N09-N16
      3 1581 1174 (N00W37 L=159)   227     85        1 (1548/1152-1599/1208)  14.0 W34-W38 N00-N01
      4 1420  985 (S08W23 L=145)   595    170        1 (1362/ 901-1459/1074)  45.0 W17-W26 S11-S03
      5 1708 1037 (S10W44 L=166)   415    135        2 (1653/1004-1748/1060) -20.6 W39-W47 S12-S08
      6 1751  949 (S16W47 L=169)   744    139        2 (1721/ 874-1772/1021)  53.1 W43-W49 S20-S12
      7  298  981 (N11E49 L= 73)   988    280        2 ( 200/ 933- 404/1019) -22.4 E57-E40 N06-N14
      8 1162  969 (S04W07 L=129)   357     91        1 (1113/ 936-1207/ 994) -31.0 W04-W09 S07-S02
      9 1580  758 (S25W31 L=153)  1212    372        0 (1415/ 700-1713/ 828)   5.0 W19-W42 S26-S22
     10 1671  675 (S32W39 L=161)  1743    323        4 (1558/ 583-1795/ 805)  25.5 W29-W50 S36-S25
     11 1061  734 (S16E03 L=119)   251    112        1 (1009/ 696-1105/ 753) -21.8 E05-W00 S18-S14
     12 1397  535 (S36W16 L=138)   238     65        0 (1370/ 509-1428/ 559)  36.9 W14-W18 S37-S34
  

    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