Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  974 1740 (N39E18 L=136)   446     93        1 ( 923/1722-1025/1757)  26.6 E21-E14 N37-N40
      2  675 1619 (N25E34 L=120)   430     96        0 ( 646/1584- 703/1659) -33.7 E36-E30 N23-N27
      3  431 1599 (N21E52 L=102)  2131    200        3 ( 360/1547- 511/1650)  -7.6 E58-E43 N17-N25
      4  480 1453 (N12E42 L=112)   930    203        1 ( 462/1355- 512/1547) -67.4 E44-E39 N06-N18
      5  283 1480 (N13E61 L= 93)  1198    253        5 ( 231/1435- 328/1529)  32.0 E67-E54 N09-N15
      6  356 1429 (N10E52 L=102)  1186    182        3 ( 313/1361- 410/1485)  54.5 E57-E47 N05-N13
      7  492 1290 (N02E37 L=117)   233     60        1 ( 469/1265- 515/1316) -26.6 E39-E35 N00-N03
      8 1567 1226 (N15W30 L=184)   229     74        0 (1541/1183-1581/1258) -53.1 W28-W31 N12-N16
      9 1061 1214 (N05E01 L=153)   312    131        0 (1008/1155-1100/1239) -36.9 E03-W02 N02-N06
     10  256 1009 (S16E52 L=102)  1103    341        1 ( 212/ 883- 351/1130)  58.0 E56-E45 S24-S08
     11 1738 1033 (N07W46 L=200)   342    117        1 (1694/ 993-1765/1059)   0.0 W42-W48 N04-N08
     12 1212  944 (S08W12 L=166)   341    144        0 (1135/ 910-1272/ 962)   6.3 W07-W16 S09-S06
     13 1358  769 (S16W25 L=179)   536    168        0 (1275/ 748-1439/ 801)  26.6 W19-W29 S17-S12
     14  483  710 (S33E33 L=121)   307     81        0 ( 450/ 685- 519/ 729)  33.7 E35-E29 S34-S30
     15  325  571 (S44E50 L=104)   832    161        1 ( 293/ 526- 361/ 614)  33.7 E54-E44 S46-S40
     16  481  493 (S48E36 L=118)  1014    215        3 ( 398/ 473- 561/ 523)   3.6 E44-E27 S49-S45
     17  862  203 (S67W07 L=161)   620     68        1 ( 837/ 179- 887/ 231) -35.5 W03-W11 S69-S63
  

    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