Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  438 1661 (N58E58 L= 30)   571     71        0 ( 413/1635- 471/1679)  14.0 E64-E51 N56-N59
      2 1554 1402 (N20W43 L=131)   742    203        0 (1507/1350-1617/1437) -18.4 W39-W48 N16-N21
      3 1268 1205 (N13W18 L=106)   235     53        0 (1232/1186-1295/1227)   0.0 W15-W19 N12-N14
      4  243 1062 (N17E57 L= 30)   450     41        3 ( 219/1039- 268/1085) -33.7 E60-E54 N16-N18
      5 1834  985 (S11W58 L=146)   713     80        2 (1798/ 961-1875/1003) -18.4 W54-W63 S12-S09
      6 1203  582 (S24W04 L= 92)  1460    593        1 ( 995/ 417-1385/ 776)  31.0 E12-W18 S32-S14
      7  995  562 (S22E09 L= 79)   874    252        0 ( 909/ 487-1070/ 657)  34.7 E16-E02 S25-S16
      8  523  545 (S17E42 L= 45)   717    173        0 ( 491/ 476- 559/ 603)  53.1 E45-E37 S21-S13
      9 1468  544 (S32W23 L=111)   391    118        0 (1425/ 494-1508/ 592)  39.8 W19-W26 S35-S28
     10 1524  490 (S38W29 L=117)   402     89        0 (1476/ 468-1562/ 514)   0.0 W24-W32 S38-S36
     11  629  421 (S27E39 L= 48)   619     74        1 ( 606/ 380- 657/ 461)  66.8 E41-E36 S30-S23
  

    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