Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  584 1711 (N30E48 L= 33)   351     39        0 ( 559/1690- 604/1730)  -9.5 E50-E44 N28-N31
      2  468 1561 (N18E49 L= 32)  1571    268        5 ( 422/1490- 546/1661) -30.1 E58-E39 N13-N25
      3 1120 1551 (N27E03 L= 77)   866    277        1 (1014/1465-1224/1618) -25.0 E11-W04 N23-N30
      4 1670 1345 (N25W38 L=119)   329     81        1 (1650/1299-1688/1384) -68.2 W36-W39 N22-N27
      5  503 1291 (N02E37 L= 43)   253     57        0 ( 472/1267- 530/1315) -21.8 E39-E34 N00-N03
      6  492 1123 (S08E35 L= 46)   804    287        1 ( 427/1048- 567/1233)  61.7 E38-E30 S13-N00
      7 1812 1135 (N16W52 L=132)   364     65        0 (1794/1107-1830/1164) -56.3 W50-W54 N13-N17
      8  698 1093 (S07E21 L= 60)   221     66        0 ( 680/1056- 716/1131)  90.0 E21-E19 S09-S04
      9 1679 1049 (N07W41 L=121)   268     56        0 (1647/1030-1710/1064)   0.0 W38-W43 N06-N08
     10  255  959 (S20E54 L= 27)   361     69        0 ( 240/ 923- 268/ 997) -53.1 E55-E52 S22-S18
     11  671  933 (S17E20 L= 60)   205     23        0 ( 646/ 922- 697/ 947)  18.4 E21-E18 S17-S16
     12  782  886 (S18E12 L= 68)   237    146        0 ( 737/ 835- 828/ 920)  50.2 E14-E09 S21-S15
     13  542  780 (S28E28 L= 52)   578    162        0 ( 498/ 716- 563/ 847)  66.0 E31-E26 S32-S23
     14  713  666 (S33E15 L= 66)  1222    224        0 ( 676/ 556- 759/ 779)  82.4 E17-E12 S40-S25
     15  391  738 (S33E41 L= 39)   255     64        0 ( 368/ 713- 422/ 769) -38.7 E43-E38 S34-S30
     16 1619  654 (S16W47 L=127)   260     45        0 (1595/ 637-1639/ 674)  33.7 W45-W48 S17-S15
     17  882  621 (S33E02 L= 78)   501    147        1 ( 804/ 605- 954/ 635)  12.5 E07-W02 S34-S31
     18  658  536 (S43E18 L= 62)   531     75        1 ( 612/ 507- 694/ 566)  23.2 E22-E15 S44-S40
     19  959  508 (S39W06 L= 86)   215     45        0 ( 929/ 484- 986/ 528)  45.0 W03-W07 S40-S37
  

    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