Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  398 1625 (N45E56 L=312)   613    165        0 ( 364/1600- 428/1640)   8.1 E59-E51 N42-N45
      2  494 1578 (N41E43 L=324)   761    130        1 ( 434/1552- 541/1604)   6.3 E48-E38 N39-N43
      3  320 1479 (N34E57 L=311)  1697    374        5 ( 274/1398- 367/1553) -34.7 E64-E49 N28-N38
      4 1445 1505 (N33W34 L= 41)   333     96        0 (1400/1485-1490/1526) -18.4 W30-W36 N31-N34
      5  784 1470 (N33E16 L=351)   243     48        1 ( 769/1447- 805/1497)  90.0 E17-E14 N31-N35
      6  860 1457 (N32E10 L=357)   366     99        2 ( 814/1443- 917/1474)  -8.1 E13-E06 N31-N33
      7 1188 1393 (N27W12 L= 20)   778    265        2 (1067/1344-1321/1419)   6.3 W03-W21 N24-N28
      8  672 1266 (N20E23 L=345)  2311    578        6 ( 471/1063- 815/1409) -38.9 E39-E12 N07-N29
      9  234 1189 (N15E59 L=309)  1623    332        4 ( 201/1068- 272/1308) -57.3 E63-E53 N07-N21
     10 1493 1211 (N14W32 L= 39)   411    134        0 (1433/1176-1548/1256) -32.0 W27-W35 N11-N16
     11 1158  851 (S07W08 L= 15)   245     84        0 (1115/ 826-1192/ 880) -36.9 W05-W09 S08-S04
     12  609  619 (S21E30 L=338)  1575    546        3 ( 436/ 489- 762/ 780)  31.0 E47-E17 S30-S10
     13 1351  740 (S15W20 L= 28)   239     46        2 (1323/ 722-1377/ 762)  18.4 W18-W21 S15-S13
     14 1731  631 (S25W52 L= 60)   425     77        1 (1705/ 614-1758/ 656)   0.0 W49-W55 S26-S23
     15  918  534 (S27E09 L=359)   691    145        2 ( 857/ 501- 999/ 572)  16.7 E13-E03 S29-S24
     16  816  470 (S32E17 L=350)   280     53        1 ( 799/ 441- 837/ 493)  53.1 E18-E15 S33-S29
  

    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