Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  396 1426 (N30E43 L= 16)  4016    490       10 ( 265/1241- 514/1617) -30.3 E78-E30 N16-N45
      2 1565 1488 (N22W43 L=102)   747    183        1 (1472/1459-1642/1510) -13.0 W35-W49 N19-N24
      3  754 1441 (N27E14 L= 45)   491    186        2 ( 697/1369- 818/1509) -54.0 E18-E09 N21-N32
      4  160 1176 (N16E65 L=354)  4553    408        8 (  85/ 995- 243/1292)  64.8 E78-E55 N05-N22
      5  225  747 (S10E62 L=357)   516     99        0 ( 206/ 702- 245/ 796)  39.8 E65-E58 S12-S07
      6  437  693 (S15E43 L= 16)   929    241        3 ( 368/ 594- 508/ 785)  33.7 E51-E35 S20-S10
      7 1293  747 (S20W14 L= 73)   222    107        1 (1251/ 716-1343/ 763)  18.4 W11-W17 S21-S19
      8 1103  669 (S23W02 L= 60)   790    241        3 (1002/ 635-1198/ 735) -32.5 E04-W07 S26-S18
      9  259  557 (S22E68 L=351)  1871    167        1 ( 218/ 468- 299/ 624)  63.4 E73-E63 S27-S17
     10 1478  465 (S41W30 L= 89)   867    110        2 (1433/ 431-1522/ 505)   0.0 W26-W34 S43-S38
     11 1294  416 (S43W15 L= 73)   242     43        1 (1278/ 397-1321/ 441) -36.9 W13-W17 S44-S40
  

    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