Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  517 1643 (N19E52 L=183)   892    136        1 ( 488/1576- 540/1702) -31.6 E59-E46 N14-N22
      2  402 1562 (N12E57 L=178)  2134    288        6 ( 363/1457- 436/1682) -47.0 E68-E49 N05-N20
      3  261 1499 (N06E70 L=164)  2638    121        1 ( 218/1448- 290/1565) -38.7 E77-E64 N03-N11
      4 1255 1326 (N16W06 L=241)   367    105        2 (1216/1284-1297/1364)  90.0 W04-W08 N12-N17
      5  277 1283 (S07E57 L=178)   880    145        2 ( 210/1255- 340/1304)   0.0 E65-E49 S07-S05
      6  241 1220 (S11E59 L=176)   448     52        1 ( 210/1204- 272/1244)   0.0 E62-E55 S12-S09
      7  363 1100 (S16E45 L=190)   252     69        1 ( 341/1072- 392/1131)  45.0 E46-E42 S18-S14
      8  343 1007 (S22E45 L=190)   843    253        1 ( 258/ 934- 428/1088) -19.4 E54-E36 S25-S18
      9  474 1061 (S17E35 L=200)   233     79        0 ( 432/1043- 506/1073)   9.5 E38-E32 S17-S15
     10  222  996 (S25E58 L=177)  1006    139        0 ( 176/ 969- 273/1032)  20.0 E64-E52 S26-S21
  

    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