Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  487 1666 (N18E55 L=277)   645     93        2 ( 460/1631- 516/1709) -71.6 E57-E53 N15-N21
      2 1146 1258 (N09W01 L=333)   487    122        0 (1096/1191-1189/1307) -29.7 E03-W04 N06-N10
      3  493 1210 (S09E36 L=297)   361     89        0 ( 458/1169- 518/1240) -36.9 E39-E34 S11-S07
      4  411 1197 (S12E42 L=291)   339     61        1 ( 385/1176- 444/1217)  45.0 E44-E39 S13-S09
      5 1539 1097 (N11W28 L=  0)   453    121        0 (1522/1030-1555/1161) -60.3 W25-W30 N07-N14
      6 1389  991 (N01W21 L=354)   745    253        0 (1308/ 876-1431/1095) -28.3 W13-W26 S04-N04
      7  689  814 (S27E14 L=318)   739    285        1 ( 562/ 786- 827/ 840)  12.5 E23-E05 S30-S24
      8  403  759 (S38E35 L=297)   254     78        0 ( 372/ 749- 440/ 775)   9.5 E38-E32 S38-S36
      9 1619  736 (S05W44 L= 17)   304     92        0 (1582/ 700-1640/ 762)  -9.5 W40-W46 S06-S04
     10  881  554 (S38W05 L=338)   311     72        0 ( 842/ 529- 918/ 584)  59.0 W03-W07 S40-S34
     11 1545  387 (S25W57 L= 30)   484    113        0 (1514/ 372-1582/ 409)   8.1 W53-W60 S26-S24
  

    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