Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1058 1833 (N48E18 L=195)  1140    235        3 (1020/1768-1097/1898) -71.6 E21-E14 N42-N57
      2  343 1507 (N13E57 L=156)   412     74        0 ( 315/1471- 369/1536)  63.4 E58-E54 N10-N14
      3 1428 1416 (N24W18 L=231)   615    147        1 (1393/1351-1463/1470)  71.6 W16-W20 N19-N28
      4  347 1384 (N05E52 L=161)   679    132        1 ( 321/1317- 380/1456)  90.0 E53-E50 N00-N09
      5  665 1395 (N09E28 L=185)   227     78        0 ( 637/1376- 708/1415)   0.0 E30-E25 N07-N10
      6 1727 1294 (N24W42 L=255)   992    146        1 (1682/1223-1760/1364) -63.4 W37-W45 N20-N28
      7  713 1228 (N00E22 L=191)   268     94        1 ( 668/1205- 764/1248)  18.4 E24-E18 S01-N01
      8  778  979 (S13E14 L=200)   995    280        4 ( 670/ 839- 865/1131)  70.7 E18-E10 S23-S03
      9 1607 1054 (N06W35 L=248)   374    116        0 (1560/1019-1650/1108) -20.6 W30-W38 N04-N08
     10  296  913 (S24E48 L=165)   749    279        0 ( 219/ 838- 356/1012)  42.5 E55-E43 S28-S17
     11  411  776 (S31E38 L=175)  1302    420        0 ( 286/ 693- 538/ 890)  35.5 E49-E28 S37-S22
     12 1590  656 (S16W44 L=257)  4617   1018        6 (1145/ 552-1870/ 738)  20.4 W17-W76 S31-S05
  

    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