Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  855 1766 (N33E37 L=264)   551     75        1 ( 830/1729- 888/1802) -66.8 E38-E34 N30-N37
      2 1214 1601 (N32E06 L=295)  1145    312        1 (1067/1554-1306/1669)  49.4 E14-W00 N26-N40
      3 1177 1380 (N17E01 L=300)   275     72        1 (1142/1348-1204/1411)  78.7 E02-W00 N14-N19
      4 1536 1265 (N21W24 L=325)   877    226        1 (1468/1212-1624/1299)  35.0 W19-W30 N15-N23
      5 1888  947 (N15W63 L=  4)  2046    232        1 (1862/ 858-1934/1018)  78.7 W59-W70 N09-N20
      6 1058  924 (S12W05 L=305)   296    119        0 (1013/ 889-1119/ 951)   0.0 W01-W08 S12-S10
      7  999  737 (S24W07 L=307)   825    306        1 ( 895/ 677-1102/ 824)  58.4 W01-W10 S29-S16
      8 1202  692 (S20W20 L=321)   382    120        1 (1146/ 657-1251/ 718)  15.9 W15-W23 S21-S18
      9 1534  624 (S14W43 L=344)   553    169        0 (1472/ 572-1606/ 661)  56.3 W40-W47 S18-S09
  

    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