Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  939 1695 (N39E23 L=308)   824    190        3 ( 866/1626-1006/1753) -20.6 E31-E15 N35-N42
      2  425 1626 (N24E61 L=270)   952     92        3 ( 401/1576- 457/1669) -33.7 E66-E55 N20-N26
      3 1370 1525 (N37W13 L=344)   369     76        0 (1329/1505-1410/1543)   0.0 W10-W16 N35-N37
      4  988 1344 (N17E08 L=323)   211    101        0 ( 952/1316-1016/1392)  71.6 E10-E07 N14-N20
      5 1450 1310 (N24W22 L=353)   324    113        0 (1424/1261-1481/1358) -45.0 W19-W24 N21-N26
      6  146 1120 (S11E69 L=262)  1382    128        2 ( 117/1076- 166/1187)  66.8 E74-E65 S14-S07
      7  180 1024 (S17E63 L=268)   982    178        1 ( 138/ 950- 217/1096)  42.3 E69-E58 S21-S11
      8 1467  926 (N01W29 L=360)   328     86        0 (1425/ 902-1503/ 946)  -9.5 W25-W31 N00-N02
      9 1864  711 (S03W72 L= 43)  1818    134        1 (1820/ 655-1895/ 765)  66.8 W67-W77 S07-N00
     10  996  677 (S23W05 L=336)  1317    429        1 ( 833/ 581-1195/ 742)  43.4 E03-W16 S31-S14
  

    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