Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  497 1487 (N33E35 L= 13)  4095    570       11 ( 308/1285- 669/1636) -29.1 E72-E20 N18-N47
      2 1696 1499 (N22W56 L=103)   970    141        1 (1628/1456-1748/1523) -26.6 W49-W61 N19-N24
      3  249 1267 (N21E55 L=352)   965    211        1 ( 224/1215- 280/1358) -71.6 E58-E51 N17-N26
      4  311 1249 (N19E48 L=359)   323     51        1 ( 294/1221- 336/1273)  45.0 E49-E45 N17-N20
      5  193 1172 (N15E60 L=347)  1529    163        3 ( 154/1116- 226/1238) -58.0 E66-E56 N11-N19
      6  213 1068 (N09E57 L=350)  1164    180        6 ( 147/1025- 279/1129) -26.6 E66-E50 N05-N13
      7 1701  995 (S08W45 L= 92)   957    177        0 (1663/ 911-1752/1046)  45.0 W41-W49 S13-S05
      8 1449  762 (S21W25 L= 72)   382    141        0 (1403/ 725-1521/ 789)  20.6 W22-W30 S22-S19
      9  511  569 (S24E41 L=  6)  1833    453        0 ( 448/ 400- 636/ 763)  35.0 E56-E26 S35-S13
     10 1272  694 (S23W13 L= 60)   674    247        1 (1181/ 660-1375/ 754) -24.8 W07-W20 S26-S18
     11  347  602 (S20E54 L=353)  1150    180        2 ( 311/ 508- 377/ 676)  74.7 E57-E50 S26-S15
     12 1574  469 (S42W41 L= 88)   328     45        0 (1550/ 453-1602/ 489) -15.9 W37-W44 S43-S40
  

    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