Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1162 1859 (N51W50 L=122)   406     53        1 (1144/1830-1181/1888)  14.0 W44-W56 N49-N53
      2  480 1668 (N62E33 L= 39)   581     31        0 ( 457/1640- 504/1692)  31.0 E36-E28 N59-N64
      3  371 1653 (N64E58 L= 14)   644     51        1 ( 345/1634- 407/1667)  -8.4 E73-E46 N62-N66
      4 1278 1520 (N25W30 L=102)   258     83        1 (1250/1485-1310/1554) -90.0 W29-W31 N22-N28
      5  852 1508 (N37W03 L= 74)   526    127        0 ( 793/1478- 917/1531) -50.2 E01-W06 N33-N39
      6  461 1474 (N46E30 L= 41)   978    153        2 ( 389/1445- 538/1498)  -3.6 E38-E22 N43-N47
      7  426 1190 (N28E35 L= 37)   690    203        1 ( 400/1138- 490/1273) -68.2 E37-E30 N23-N33
      8  180 1118 (N28E62 L= 10)   667    118        0 ( 159/1065- 198/1167) -50.2 E65-E59 N25-N31
      9  777 1126 (N16E12 L= 60)   232     99        1 ( 764/1071- 797/1162)  45.0 E13-E09 N13-N17
     10 1519 1035 (S08W29 L=101)   493    128        2 (1488/ 980-1543/1096)  45.0 W25-W31 S10-S04
     11  587 1030 (N15E26 L= 46)   540    139        3 ( 547/ 970- 628/1074) -81.9 E28-E22 N11-N18
     12 1486  736 (S24W21 L= 92)   839    212        2 (1454/ 668-1531/ 847)  85.2 W18-W23 S29-S17
     13 1685  683 (S34W36 L=108)   649     87        3 (1662/ 638-1715/ 730) -60.3 W34-W39 S37-S30
     14 1282  627 (S25W05 L= 77)   299     84        0 (1237/ 616-1328/ 638) -31.0 W02-W08 S26-S23
     15  619  318 (S28E51 L= 20)  2269    288        1 ( 553/ 203- 725/ 418)  16.3 E77-E36 S36-S23
     16 1451  340 (S53W12 L= 83)   521     52        0 (1423/ 322-1486/ 364) -39.8 W09-W15 S55-S49
     17 1311  229 (S58E06 L= 65)   648     46        1 (1285/ 205-1335/ 256) -90.0 E10-E03 S62-S54
     18 1219  176 (S60E20 L= 52)   569     54        1 (1197/ 154-1248/ 196)  12.5 E23-E14 S62-S56
  

    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