Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1224 1498 (N26W08 L=319)   533    169        1 (1169/1421-1253/1563) -48.8 W03-W10 N21-N29
      2  488 1335 (N09E37 L=273)  1892    380        3 ( 434/1175- 566/1490) -69.8 E45-E31 S01-N18
      3 1694 1178 (N12W42 L=353)  1976    457        5 (1591/1068-1832/1259) -21.8 W33-W55 N07-N16
      4 1169 1183 (N05W07 L=318)   318    118        1 (1133/1145-1214/1207)   0.0 W04-W09 N03-N06
      5 1272 1117 (N03W14 L=325)   303    147        1 (1230/1051-1297/1160) -59.0 W10-W15 S01-N05
      6  646  976 (S12E23 L=288)   293     63        1 ( 609/ 954- 689/ 996)   9.5 E25-E19 S13-S10
      7 1575  834 (S10W37 L=348)   477    138        0 (1525/ 800-1620/ 891) -24.0 W32-W41 S11-S07
      8  178  828 (S23E61 L=250)   515    137        3 ( 154/ 800- 200/ 856) -71.6 E64-E58 S24-S21
      9 1206  655 (S26W16 L=327)  1646    388        1 (1031/ 570-1396/ 792)  39.0 W05-W26 S32-S15
     10  900  537 (S37E03 L=308)   404     92        2 ( 852/ 521- 940/ 564)  26.6 E06-W00 S38-S34
  

    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