Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1262 1590 (N39W04 L=149)   269     49        0 (1238/1572-1296/1606)   0.0 W02-W07 N37-N40
      2  926 1472 (N24E16 L=129)   634    238        1 ( 843/1384-1009/1552) -28.1 E23-E08 N19-N27
      3 1376 1532 (N38W14 L=159)   280     51        0 (1357/1510-1403/1549)  36.9 W12-W16 N35-N39
      4 1492 1486 (N37W24 L=169)   260     63        0 (1464/1470-1527/1501)  31.0 W21-W26 N35-N38
      5 1670 1382 (N34W39 L=184)  2073    363        2 (1539/1311-1793/1468)  -9.9 W27-W51 N30-N37
      6 1392 1168 (N14W20 L=164)  1207    424        4 (1336/1050-1479/1261) -39.3 W15-W26 N07-N20
      7 1133  990 (S01W07 L=152)   368    131        1 (1056/ 959-1174/1033)  39.8 W02-W09 S04-N01
      8  483  706 (S31E31 L=114)   504    168        0 ( 438/ 646- 513/ 768)  63.4 E34-E28 S36-S26
      9 1137  700 (S18W13 L=158)   245     40        0 (1117/ 675-1159/ 724) -26.6 W11-W15 S19-S16
     10 1602  666 (S11W46 L=190)   833    191        1 (1536/ 636-1668/ 699)   0.0 W39-W51 S12-S09
     11 1487  384 (S31W50 L=195)   631    131        1 (1452/ 365-1536/ 407)   0.0 W45-W54 S32-S29
     12 1236  255 (S47W40 L=185)   355     53        0 (1212/ 238-1254/ 280) -76.0 W38-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