Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1488 1610 (N24W51 L=179)   521    103        1 (1438/1580-1536/1633) -51.3 W47-W52 N20-N26
      2 1364 1471 (N21W34 L=163)  1240    493        0 (1146/1424-1555/1533) -36.5 W23-W46 N13-N30
      3 1802 1333 (S01W62 L=190)   526     74        1 (1769/1304-1833/1361) -38.7 W59-W64 S02-N02
      4 1142 1244 (N14W13 L=141)   242    117        0 (1089/1229-1201/1253) -26.6 W09-W15 N12-N15
      5  811 1189 (N20E08 L=121)   421    148        2 ( 743/1167- 878/1209) -26.6 E12-E04 N17-N22
      6 1902 1053 (S19W64 L=193)  1302    161        3 (1877/ 980-1926/1129)  77.5 W59-W68 S23-S14
      7  315 1090 (N26E46 L= 82)   387     61        3 ( 283/1065- 351/1110) -15.9 E49-E42 N24-N27
      8  187 1069 (N28E62 L= 67)   393     52        3 ( 167/1046- 206/1093) -31.0 E64-E58 N25-N29
      9  525 1005 (N17E30 L= 98)   221    101        0 ( 475/ 997- 576/1015) -14.0 E34-E26 N15-N17
     10 1298  638 (S24W05 L=134)   971    270        1 (1260/ 521-1316/ 769)  57.3 W00-W09 S30-S15
     11  506  509 (S11E48 L= 81)   291     49        0 ( 489/ 479- 526/ 538)  63.4 E49-E46 S13-S09
     12 1249  472 (S33E02 L=126)   267     43        1 (1225/ 453-1279/ 491) -18.4 E03-W00 S34-S31
  

    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