Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  840 1796 (N35E40 L=299)   617     88        0 ( 817/1755- 866/1838) -58.0 E43-E36 N31-N39
      2  640 1737 (N26E50 L=289)   889    190        0 ( 556/1698- 727/1766)  37.9 E55-E43 N21-N30
      3 1037 1518 (N21E13 L=326)   726    210        0 ( 953/1463-1119/1543)   5.2 E17-E06 N18-N23
      4  813 1423 (N10E23 L=316)   623    234        0 ( 777/1329- 845/1528) -42.3 E28-E17 N05-N15
      5  639 1357 (N02E31 L=308)   337     69        0 ( 621/1328- 670/1392) -26.6 E32-E28 N00-N04
      6 1735 1313 (N31W40 L= 19)  1387    209        2 (1666/1233-1796/1373) -17.1 W32-W46 N27-N34
      7 1632 1095 (N14W34 L= 13)  1113    349        1 (1543/ 959-1707/1230) -29.1 W25-W43 N08-N19
      8 1459 1135 (N11W22 L=  0)   742    193        0 (1435/1052-1492/1220) -73.3 W20-W24 N06-N16
      9  493 1160 (S12E35 L=303)   413    137        0 ( 441/1144- 565/1183)  18.4 E39-E30 S13-S10
     10  231 1111 (S20E58 L=281)   776    134        1 ( 192/1055- 262/1173)  58.0 E61-E54 S23-S15
     11  520  989 (S21E30 L=309)   285    110        0 ( 509/ 941- 546/1035)  78.7 E30-E28 S23-S18
     12  668  952 (S20E19 L=320)   470    146        2 ( 599/ 933- 731/ 983)  32.0 E23-E15 S21-S16
     13  446  893 (S29E34 L=305)   333    117        0 ( 398/ 876- 504/ 917)  24.0 E38-E29 S30-S25
     14  924  842 (S20E01 L=338)   190     71        0 ( 888/ 818- 956/ 857) -14.0 E03-W01 S20-S18
     15  665  733 (S33E15 L=324)   462     93        2 ( 621/ 701- 710/ 763)  45.0 E17-E11 S35-S30
     16  586  697 (S37E20 L=319)   386     81        0 ( 556/ 671- 626/ 723)  -9.5 E22-E16 S38-S35
     17 1367  466 (S28W41 L= 19)  1125    180        2 (1281/ 427-1444/ 513)  59.0 W36-W44 S32-S22
     18 1147  318 (S44W37 L= 16)   493     87        0 (1103/ 303-1193/ 332)  15.9 W32-W39 S45-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