Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1613 1713 (N25W75 L= 42)   903     67        2 (1581/1687-1638/1731)  -9.5 W70-W83 N23-N26
      2 1437 1596 (N26W46 L= 13)   449     96        0 (1382/1571-1485/1618) -45.0 W43-W49 N23-N28
      3 1376 1491 (N22W36 L=  3)   589    168        1 (1296/1445-1441/1522) -53.1 W32-W38 N17-N26
      4 1153 1323 (N20W16 L=342)  1117    425        9 (1031/1187-1299/1476) -77.7 W12-W19 N08-N31
      5  314 1436 (N50E50 L=276)   522     54        0 ( 296/1409- 340/1461)  45.0 E53-E46 N47-N51
      6  553 1320 (N36E23 L=303)   963    216        2 ( 470/1287- 653/1357) -29.7 E29-E15 N31-N40
      7 1639  984 (S13W36 L=  2)  1156    292        0 (1578/ 858-1691/1114)  38.2 W28-W43 S18-S07
      8  894 1017 (N09E07 L=319)   578    188        0 ( 829/ 948- 968/1082) -70.0 E09-E04 N03-N14
      9  281  967 (N21E52 L=275)  1309    311        3 ( 233/ 879- 322/1080) -90.0 E56-E48 N14-N27
     10  598  951 (N13E27 L=300)   325     68        0 ( 570/ 920- 637/ 978) -53.1 E28-E24 N10-N14
     11 1484  853 (S16W22 L=349)   255    108        0 (1433/ 835-1537/ 874) -33.7 W19-W25 S18-S13
     12 1447  735 (S22W17 L=344)   361    149        0 (1384/ 710-1508/ 777)  -5.7 W12-W22 S22-S19
     13 1645  717 (S29W31 L=358)   836    118        3 (1602/ 676-1692/ 756)  14.0 W27-W36 S31-S26
     14  270  732 (N07E59 L=267)   448     80        1 ( 235/ 718- 311/ 747)   0.0 E64-E54 N06-N08
     15  851  339 (S30E33 L=293)   676    119        2 ( 808/ 300- 902/ 376)  16.7 E37-E27 S31-S28
     16 1102  309 (S40E19 L=308)   636    123        1 (1070/ 256-1120/ 372)  48.4 E23-E15 S44-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