Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  697 1550 (N42E10 L=283)   323     76        0 ( 670/1515- 733/1576)  53.1 E12-E06 N39-N43
      2 1767 1530 (N14W68 L=  1)  1357    279        4 (1742/1469-1799/1572)  24.4 W62-W74 N10-N16
      3  821 1439 (N31E03 L=290)   314    105        0 ( 764/1425- 875/1451) -15.9 E07-W00 N29-N32
      4 1279 1249 (N10W20 L=313)   841    132        4 (1189/1218-1388/1277) -28.6 W14-W25 N05-N12
      5  871 1193 (N16E05 L=288)  1353    329        2 ( 703/1148-1031/1226) -11.3 E16-W04 N11-N18
      6  952 1095 (N08E03 L=291)   180     87        0 ( 919/1053- 986/1132)  31.0 E05-W00 N06-N09
      7 1744  963 (S17W45 L=338)   644    123        0 (1714/ 885-1769/1021)  45.0 W41-W47 S21-S13
      8  133  993 (N18E65 L=228)   936    226        0 (  89/ 968- 185/1014) -14.0 E71-E59 N16-N19
      9  380  985 (N13E41 L=252)   319     72        2 ( 348/ 954- 404/1009) -36.9 E43-E39 N11-N15
     10  280  840 (N07E53 L=240)   598     87        2 ( 253/ 811- 315/ 868)   0.0 E56-E49 N05-N08
     11 1352  797 (S17W14 L=308)   267    121        0 (1342/ 733-1376/ 849)  66.8 W12-W16 S21-S14
     12 1211  651 (S23W03 L=296)  1320    290        3 (1130/ 505-1291/ 759)  54.5 E05-W08 S30-S15
     13  560  587 (S13E39 L=254)  1090    322        5 ( 463/ 476- 651/ 667) -74.1 E44-E35 S20-S06
     14  411  562 (S11E51 L=242)  1985    382        6 ( 337/ 499- 474/ 627)  12.1 E59-E45 S15-S07
     15  308  532 (S11E63 L=230)  2043    281        7 ( 262/ 468- 362/ 599) -63.4 E67-E58 S15-S07
     16  674  550 (S17E33 L=261)   385     90        1 ( 643/ 505- 700/ 586)  51.3 E34-E29 S20-S15
     17  555  502 (S18E43 L=250)   407    108        1 ( 492/ 486- 599/ 514) -14.0 E47-E39 S19-S16
     18 1079  325 (S42E16 L=277)   554    163        2 (1050/ 281-1106/ 375)  69.4 E18-E12 S45-S37
     19 1001  200 (S51E30 L=263)   971    190        2 ( 972/ 158-1033/ 248)  52.1 E35-E24 S55-S46
  

    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