Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  660 1677 (N40E41 L=124)   615     96        1 ( 619/1646- 695/1717) -18.4 E48-E36 N37-N42
      2  442 1609 (N32E58 L=106)   613     73        0 ( 424/1572- 464/1650) -24.4 E64-E53 N29-N34
      3 1294 1581 (N41W14 L=179)   757    149        3 (1231/1543-1360/1623) -15.3 W09-W20 N38-N43
      4  659 1424 (N22E30 L=135)   566    119        3 ( 626/1368- 692/1482) -82.9 E31-E27 N18-N26
      5 1517 1206 (N17W31 L=195)   341    161        0 (1494/1161-1551/1265) -63.4 W29-W33 N14-N20
      6  304  870 (S16E50 L=114)  3963    490        6 ( 198/ 659- 409/1050)  43.0 E72-E40 S31-S03
      7 1711  727 (S11W52 L=217)   545    110        0 (1696/ 678-1725/ 784) -74.1 W50-W54 S13-S06
      8 1672  717 (S12W49 L=214)   347    134        0 (1630/ 674-1690/ 774)  74.1 W46-W50 S14-S07
      9 1265  743 (S14W18 L=183)   316     92        0 (1219/ 729-1312/ 762)   9.5 W15-W21 S14-S12
     10 1248  555 (S26W21 L=185)   337     79        1 (1205/ 533-1288/ 579)  31.0 W18-W23 S27-S24
     11 1634  501 (S26W55 L=220)  1565    164        2 (1561/ 460-1706/ 543)  22.6 W49-W62 S29-S22
     12 1277  434 (S34W27 L=191)   711    142        0 (1206/ 414-1338/ 455)  11.3 W20-W30 S35-S32
  

    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