Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  470 1551 (N29E50 L=  2)   435     68        0 ( 455/1513- 496/1575)  76.0 E52-E48 N26-N31
      2  752 1441 (N25E23 L= 29)  1297    323        2 ( 653/1367- 807/1558)  77.0 E28-E18 N20-N33
      3 1671 1351 (N28W46 L= 98)  3245    547        1 (1482/1172-1817/1503) -34.2 W30-W59 N16-N36
      4 1550 1160 (N14W34 L= 86)   562    120        1 (1514/1123-1593/1209) -45.0 W30-W37 N12-N17
      5  215  947 (S11E60 L=353)   643     62        1 ( 184/ 922- 250/ 969)   0.0 E63-E55 S12-S09
      6  615  758 (S19E25 L= 28)   516    146        1 ( 589/ 689- 654/ 811)  63.4 E26-E21 S23-S15
      7 1623  646 (S17W46 L= 99)  1098    219        0 (1537/ 576-1696/ 720)  48.0 W41-W51 S22-S11
      8  697  640 (S26E19 L= 34)   259     52        0 ( 670/ 620- 725/ 659) -26.6 E20-E16 S26-S24
      9 1439  492 (S29W36 L= 89)   443     97        1 (1391/ 471-1488/ 517)  23.2 W32-W40 S30-S26
  

    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