Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1697 1444 (N20W56 L=105)  1524    391        5 (1573/1321-1768/1583) -81.0 W51-W61 N11-N30
      2  753 1377 (N32E13 L= 35)   268     97        0 ( 706/1366- 805/1389) -15.9 E16-E09 N30-N32
      3  417 1291 (N30E40 L=  8)  1977    388        5 ( 283/1220- 563/1341)   4.8 E52-E28 N25-N32
      4 1361 1261 (N15W25 L= 73)   273     88        1 (1332/1216-1379/1307) -90.0 W23-W25 N12-N18
      5  169 1231 (N27E65 L=344)   765    230        3 ( 134/1193- 197/1265) -29.7 E69-E60 N24-N28
      6  883 1154 (N16E07 L= 41)   739    234        3 ( 773/1092- 969/1235) -47.7 E13-E02 N11-N22
      7  340  645 (S10E54 L=354)  3171    617        6 ( 273/ 479- 382/ 808)  49.2 E71-E47 S22-N00
      8  546  725 (S06E35 L= 13)   736    173        4 ( 479/ 691- 620/ 753) -24.0 E39-E30 S09-S04
      9 1452  666 (S23W23 L= 71)   266     32        2 (1430/ 645-1475/ 690)  33.7 W21-W24 S23-S20
     10  489  635 (S12E41 L=  7)   531    114        1 ( 472/ 572- 510/ 686)  81.9 E43-E39 S15-S08
     11  458  366 (S31E56 L=352)  1751    345        2 ( 397/ 330- 542/ 414)  14.7 E64-E45 S34-S27
     12 1170  336 (S43E01 L= 48)   309     48        0 (1139/ 313-1192/ 356)  36.9 E03-W01 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