Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1513 1162 (N16W28 L=284)  4886   1038       15 (1187/ 781-1747/1495) -31.6 W01-W54 S02-N31
      2  291 1399 (N08E61 L=195)   546    103        0 ( 272/1374- 321/1432) -63.4 E63-E57 N06-N10
      3  320 1014 (S14E47 L=208)   452     61        1 ( 296/ 984- 344/1045)  76.0 E49-E45 S15-S11
      4  213  955 (S19E57 L=198)   520    114        1 ( 185/ 911- 236/ 991)  56.3 E60-E54 S22-S16
      5  345  896 (S21E44 L=212)   657    217        1 ( 284/ 838- 396/ 981)  54.0 E49-E40 S25-S14
      6 1587  736 (S08W42 L=297)   780    205        3 (1515/ 713-1677/ 761)  22.6 W36-W49 S10-S05
      7  367  717 (S32E42 L=214)   383    103        1 ( 332/ 685- 411/ 747) -23.2 E44-E37 S33-S30
      8  737  710 (S25E13 L=243)   468    134        1 ( 665/ 684- 797/ 729)   0.0 E17-E08 S26-S23
      9 1261  640 (S20W22 L=278)   559    203        0 (1174/ 617-1354/ 690)   0.0 W15-W28 S20-S18
     10 1294  282 (S43W42 L=297)   383     62        1 (1267/ 262-1314/ 301)  71.6 W39-W44 S44-S41
  

    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