Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1004 1752 (N51W29 L= 10)   301     45        1 ( 980/1729-1030/1767) -11.3 W26-W32 N49-N52
      2 1532 1599 (N23W54 L= 35)  1004    190        2 (1501/1495-1553/1665)  29.7 W45-W60 N17-N26
      3 1300 1505 (N25W32 L= 13)   455    134        0 (1242/1469-1351/1529) -60.3 W29-W34 N21-N28
      4 1197 1396 (N22W21 L=  2)   662    222        2 (1097/1343-1273/1433) -45.0 W15-W24 N17-N26
      5  240 1403 (N49E63 L=278)   483     59        0 ( 215/1379- 271/1418)  11.3 E67-E56 N47-N50
      6  970 1207 (N18W02 L=343)   891    361        2 ( 886/1074-1095/1328) -78.7 E02-W05 N07-N27
      7  404 1257 (N36E37 L=304)  1101    190        3 ( 328/1208- 494/1307) -29.4 E45-E29 N31-N40
      8  572 1132 (N24E25 L=316)   313     74        0 ( 541/1100- 605/1161) -59.0 E26-E22 N21-N26
      9  171 1023 (N26E65 L=275)   743     69        3 ( 154/ 983- 186/1061)  63.4 E68-E62 N23-N27
     10  196  931 (N20E63 L=278)  2468    210        7 ( 135/ 853- 244/1036)  51.3 E72-E55 N15-N25
     11 1741  938 (S19W43 L= 24)   878    176        2 (1721/ 857-1758/1033)  68.2 W40-W46 S23-S13
     12 1459  884 (S13W21 L=  2)   951    298        0 (1399/ 769-1517/1026)  38.2 W14-W28 S18-S06
     13  697  913 (N09E22 L=319)   566    172        0 ( 637/ 854- 762/ 978) -66.0 E23-E18 N04-N13
     14  408  886 (N13E42 L=299)   295     73        0 ( 373/ 859- 430/ 905) -31.0 E45-E40 N11-N14
     15 1209  627 (S21E00 L=341)   299    117        0 (1175/ 608-1260/ 659) -45.0 E01-W03 S22-S18
     16 1458  609 (S30W15 L=356)   587    111        1 (1424/ 557-1496/ 644)  29.7 W11-W18 S32-S27
     17  974  449 (S26E20 L=321)   434    171        0 ( 927/ 392-1035/ 512)  22.6 E25-E13 S28-S23
     18  954  262 (S39E32 L=309)   690    126        0 ( 925/ 204- 977/ 325)  54.0 E37-E28 S44-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