Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  360 1498 (N33E47 L=326)   542    111        1 ( 337/1429- 388/1540)  90.0 E48-E45 N27-N35
      2 1876 1325 (N12W67 L= 81)   552     86        3 (1851/1289-1896/1359) -63.4 W64-W70 N10-N14
      3  308 1254 (N17E48 L=325)   412     73        1 ( 287/1220- 326/1286) -63.4 E50-E46 N14-N18
      4 1057 1237 (N10W03 L= 17)   250     50        0 (1037/1209-1071/1274)  76.0 W01-W04 N08-N12
      5 1454 1238 (N08W28 L= 41)   231     74        0 (1416/1223-1491/1249) -21.8 W25-W30 N07-N09
      6  586 1136 (N07E26 L=347)   948    278        1 ( 470/1035- 669/1204) -37.6 E34-E21 N00-N11
      7  279 1031 (N03E50 L=323)   967    198        4 ( 249/ 960- 319/1083)  63.4 E53-E46 N00-N06
      8  368 1026 (N02E43 L=331)   373     96        0 ( 332/ 999- 411/1051) -18.4 E45-E39 N00-N04
      9 1643  680 (S26W42 L= 55)   339     69        0 (1617/ 639-1675/ 710)  53.1 W39-W43 S28-S24
     10 1179  446 (S40W08 L= 21)   247     42        0 (1151/ 431-1208/ 460)  14.0 W05-W09 S40-S38
  

    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