Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  599 1322 (N30E25 L=208)   338     94        0 ( 552/1302- 645/1343) -23.2 E28-E21 N28-N32
      2 1194 1195 (N14W13 L=246)   506    213        0 (1087/1134-1275/1224) -38.7 W06-W17 N09-N17
      3  166 1197 (N25E64 L=168)  1045    218        1 ( 123/1177- 205/1214)  10.3 E70-E59 N24-N26
      4  422 1053 (N15E40 L=193)  1865    386        4 ( 373/ 858- 494/1204) -79.7 E43-E35 N02-N24
      5  516 1167 (N21E32 L=201)   287     60        0 ( 486/1140- 540/1193) -45.0 E33-E29 N19-N23
      6 1127  954 (N01W05 L=238)   341    104        0 (1073/ 919-1170/ 989) -38.7 W02-W07 S01-N03
      7  254  856 (N04E58 L=175)  1182    265        6 ( 202/ 780- 296/ 942)  48.4 E64-E52 N00-N09
      8  995  532 (S24E10 L=223)   470    134        3 ( 970/ 474-1011/ 592)  77.5 E11-E07 S28-S19
      9 1414  255 (S58W24 L=256)   611     83        2 (1374/ 238-1452/ 274)   0.0 W18-W29 S60-S55
  

    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