Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1097 1602 (N35W20 L=130)   793    209        1 ( 986/1576-1184/1623) -11.3 W11-W26 N33-N37
      2 1619 1521 (N17W54 L=163)   469     63        0 (1597/1487-1644/1552)  23.2 W50-W57 N15-N18
      3 1712 1520 (N15W63 L=173)   487     30        1 (1687/1502-1737/1543) -21.8 W59-W65 N13-N16
      4  415 1309 (N32E36 L= 73)  2064    312        5 ( 298/1212- 533/1378) -30.6 E48-E26 N24-N38
      5  897 1253 (N18E03 L=107)  1266    441        2 ( 707/1203-1070/1368)   2.4 E14-W10 N15-N21
      6  225 1191 (N29E56 L= 54)   560     58        2 ( 205/1160- 247/1219)  21.8 E58-E53 N26-N30
      7 1815  899 (S23W53 L=162)  2373    446        4 (1734/ 733-1885/1098)  48.2 W43-W63 S31-S12
      8 1527  744 (S25W26 L=135)   651    163        1 (1486/ 675-1567/ 822)  48.8 W22-W29 S28-S20
      9  662  631 (S12E31 L= 79)   487    133        1 ( 644/ 561- 679/ 692)  63.4 E33-E29 S16-S08
     10 1206  613 (S25W02 L=111)   307     31        0 (1182/ 590-1229/ 639)  33.7 W00-W03 S26-S23
     11  437  516 (S15E52 L= 57)   822    132        1 ( 407/ 461- 483/ 566)  76.0 E54-E48 S19-S11
     12 1154  481 (S33E05 L=105)   279     44        1 (1131/ 453-1173/ 503) -90.0 E05-E03 S34-S30
  

    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