Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  675 1679 (N33E33 L= 94)   371     50        2 ( 653/1653- 698/1707)  90.0 E35-E31 N30-N35
      2 1149 1550 (N28W03 L=130)  1207    436        4 (1038/1480-1336/1678) -25.5 E06-W15 N24-N36
      3 1620 1270 (N16W36 L=163)   335     53        3 (1596/1239-1641/1296)  56.3 W34-W37 N14-N17
      4 1726 1228 (N15W45 L=172)   297     91        0 (1706/1185-1757/1270) -51.3 W43-W47 N12-N17
      5 1549 1126 (N06W31 L=159)  1326    395        0 (1410/ 977-1619/1268) -43.0 W21-W37 S01-N13
      6  352 1149 (S03E45 L= 83)   256     69        1 ( 326/1119- 383/1186) -31.0 E47-E42 S04-N00
      7 1025 1068 (S03E00 L=127)   432    170        0 ( 946/1016-1084/1094) -14.0 E05-W03 S05-S01
      8  166  892 (S19E64 L= 64)  4794    370        5 (  95/ 761- 215/1068)  52.1 E83-E57 S26-S08
      9  104 1004 (S12E72 L= 55)   776     84        1 (  71/ 980- 133/1031)  -8.7 E80-E67 S13-S09
     10  358  892 (S18E44 L= 84)  1081    313        2 ( 293/ 826- 418/ 967) -45.0 E49-E38 S22-S13
     11 1803  866 (S05W55 L=183)   894    110        5 (1757/ 825-1855/ 905)  -5.2 W50-W61 S07-S03
     12 1282  752 (S19W19 L=146)   848    357        1 (1214/ 642-1384/ 884) -40.6 W13-W27 S24-S11
     13 1741  727 (S14W53 L=181)  1824    199        4 (1667/ 658-1819/ 809)  73.3 W47-W60 S18-S08
     14 1551  725 (S17W38 L=165)  1203    346        1 (1410/ 666-1671/ 767)  27.9 W29-W46 S22-S13
     15 1222  689 (S23W16 L=144)   203     59        0 (1207/ 651-1238/ 722)  90.0 W15-W17 S25-S21
     16  972  643 (S29W01 L=128)   248     67        0 ( 949/ 607- 999/ 681)  76.0 W00-W02 S31-S26
     17  322  549 (S40E56 L= 71)  1856    156        5 ( 243/ 499- 394/ 586)   7.9 E75-E46 S43-S37
     18 1636  393 (S35W63 L=191)  1250     92        4 (1595/ 360-1686/ 429)   3.8 W56-W71 S37-S33
  

    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