Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  980 1814 (N52E25 L=199)   470     59        0 ( 957/1793-1012/1832)  -8.1 E28-E21 N49-N54
      2 1045 1723 (N46E14 L=211)  1155    151        5 ( 990/1685-1107/1775)   5.2 E18-E07 N42-N50
      3 1055 1602 (N36E09 L=216)   501     81        1 (1016/1575-1095/1644)  18.4 E11-E05 N33-N39
      4  369 1542 (N20E60 L=165)  1382    231        0 ( 300/1445- 443/1618) -15.7 E80-E48 N15-N24
      5  279 1385 (N09E61 L=164)  1792    221        3 ( 234/1294- 322/1484) -80.5 E65-E57 N03-N15
      6 1345 1173 (N13W17 L=242)   787    176        1 (1278/1128-1404/1215)  45.0 W13-W20 N09-N16
      7 1297 1102 (N08W15 L=240)   947    269        3 (1221/1029-1391/1209) -22.6 W09-W21 N04-N13
      8  784 1171 (N04E17 L=208)   240     85        0 ( 748/1149- 826/1198)  -9.5 E19-E13 N03-N05
      9 1068 1117 (N06W01 L=226)   436    133        0 (1050/1052-1088/1186) -74.1 W00-W03 N01-N09
     10  404 1033 (S10E40 L=185)   724    231        0 ( 388/ 911- 427/1123) -85.2 E42-E38 S16-S04
     11 1568 1009 (N07W34 L=259)   293    134        0 (1517/ 970-1601/1035)  -8.1 W30-W37 N05-N08
     12  315  855 (S22E47 L=178)  1088    188        2 ( 223/ 812- 400/ 881)  14.7 E58-E39 S26-S19
     13 1296  739 (S13W21 L=246)   697    251        2 (1242/ 642-1363/ 849)  77.9 W19-W23 S19-S05
     14 1675  608 (S15W54 L=279)   325     59        0 (1638/ 589-1697/ 625)  -8.1 W49-W56 S15-S13
     15 1557  439 (S27W53 L=277)   835    180        2 (1469/ 398-1636/ 478)  30.3 W46-W58 S31-S23
     16  492  420 (S50E36 L=188)   570     97        1 ( 447/ 405- 539/ 439)  12.1 E44-E30 S52-S48
     17  598  389 (S51E25 L=200)   414     58        0 ( 580/ 366- 620/ 421) -68.2 E26-E22 S52-S47
     18  694  309 (S56E15 L=210)   937    111        1 ( 656/ 274- 740/ 354) -54.5 E19-E09 S58-S51
     19  970  235 (S56W17 L=241)  1080    150        2 ( 925/ 182-1026/ 304) -26.6 W08-W26 S60-S50
  

    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