Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  916 1830 (N64E15 L=  9)  1064    163        2 ( 849/1807- 994/1852)   0.0 E23-E03 N61-N66
      2  849 1610 (N43E15 L=  9)   257     63        1 ( 833/1573- 865/1644) -63.4 E16-E13 N39-N45
      3 1366 1411 (N28W25 L= 48)   303     78        0 (1331/1386-1405/1439) -31.0 W22-W27 N26-N29
      4 1030 1324 (N23W00 L= 24)   304    123        0 ( 985/1278-1077/1363) -45.0 E02-W03 N19-N25
      5  326 1288 (N19E52 L=331)   360     61        0 ( 298/1265- 360/1313)  31.0 E54-E49 N17-N20
      6 1822  787 (S13W62 L= 85)  1077    262        1 (1779/ 668-1845/ 910)  82.4 W59-W64 S20-S05
      7 1615  854 (S08W40 L= 63)   259    110        0 (1589/ 820-1638/ 878)  45.0 W37-W41 S09-S06
      8 1654  791 (S12W44 L= 67)  1851    384        3 (1532/ 707-1774/ 878)  26.6 W34-W54 S17-S06
      9 1071  699 (S17W03 L= 26)   718    193        2 ( 983/ 661-1155/ 736)   5.7 E02-W08 S19-S14
     10 1548  581 (S26W38 L= 62)  1138    209        1 (1454/ 547-1614/ 650)   4.4 W30-W44 S28-S21
     11  521  549 (S28E38 L=346)   436     69        0 ( 492/ 524- 543/ 575)  31.0 E40-E35 S29-S25
     12  529  473 (S33E40 L=344)  1840    267        4 ( 445/ 413- 625/ 569) -39.8 E44-E32 S38-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