Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1047 1649 (N39E13 L=194)   567    103        0 (1019/1606-1084/1692) -29.7 E16-E09 N36-N42
      2 1347 1473 (N33W13 L=220)   802    149        2 (1289/1409-1408/1526) -33.7 W08-W18 N29-N36
      3  376 1492 (N15E57 L=150)   372     62        2 ( 352/1464- 398/1524)  90.0 E59-E55 N13-N17
      4 1563 1182 (N19W32 L=239)   413    153        0 (1542/1113-1606/1246) -49.4 W29-W35 N15-N22
      5 1340 1110 (N10W17 L=225)   359    109        1 (1306/1082-1362/1151)  78.7 W15-W19 N07-N12
      6 1840 1025 (N15W60 L=267)   693     71        1 (1813/ 995-1866/1055) -21.8 W56-W63 N13-N16
      7  379  893 (S21E41 L=166)  2364    429        5 ( 221/ 761- 511/1029)  38.2 E60-E32 S32-S10
      8 1268  673 (S17W22 L=229)   316     41        1 (1243/ 650-1293/ 695) -14.0 W20-W24 S18-S16
      9 1739  621 (S11W61 L=269)  1535    157        2 (1660/ 592-1795/ 661)  22.6 W53-W71 S14-S08
     10  940  265 (S55W15 L=223)   837    111        1 ( 877/ 237- 986/ 297)   9.5 W08-W20 S56-S52
     11 1172  262 (S48W36 L=243)   446     75        1 (1143/ 242-1217/ 288)   0.0 W31-W41 S49-S46
     12  757  260 (S61E04 L=204)   483     80        0 ( 715/ 246- 795/ 280)   6.3 E09-W00 S63-S58
  

    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