Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  519 1511 (N31E33 L=326)   931    251        1 ( 477/1381- 568/1608)  79.4 E35-E31 N22-N38
      2 1550 1342 (N15W36 L= 36)   407    104        3 (1492/1313-1586/1372) -38.7 W32-W38 N12-N16
      3 1279 1275 (N11W17 L= 17)   200     43        0 (1263/1250-1295/1302) -90.0 W15-W17 N10-N13
      4  846 1145 (N06E10 L=350)   801    224        1 ( 760/1035- 899/1228) -51.3 E15-E07 N00-N11
      5  453 1099 (N06E36 L=324)   307     69        1 ( 420/1077- 484/1123) -11.3 E38-E33 N04-N07
      6  449 1034 (N02E36 L=323)   447    105        1 ( 416/1006- 496/1073) -38.7 E38-E33 N00-N04
      7 1781  680 (S26W54 L= 54)   644    121        0 (1747/ 622-1818/ 743)  63.4 W51-W57 S30-S22
      8 1403  507 (S37W25 L= 25)   719    170        2 (1325/ 482-1473/ 541)   0.0 W18-W30 S38-S34
      9 1230  392 (S44W12 L= 12)   544    142        0 (1157/ 368-1299/ 424)  10.3 W06-W18 S46-S42
  

    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