Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  758 1874 (N39E59 L=350)   929     75        0 ( 743/1833- 774/1913) -25.2 E70-E52 N35-N43
      2 1061 1641 (N31E18 L= 31)   307     67        1 (1036/1608-1084/1672) -80.5 E19-E15 N27-N33
      3  350 1563 (N08E65 L=344)  3678    308        3 ( 259/1468- 484/1662)  53.7 E72-E57 N01-N16
      4  688 1558 (N15E36 L= 12)   250     99        0 ( 661/1523- 727/1591) -23.2 E39-E32 N13-N16
      5  953 1250 (N04E10 L= 38)   352    140        0 ( 927/1192- 997/1313) -50.2 E13-E06 N01-N07
      6  149 1184 (S18E72 L=337)  1035     71        0 ( 123/1152- 174/1219)  33.7 E77-E67 S19-S15
      7 1588 1152 (N18W29 L= 78)   248     69        0 (1569/1118-1609/1185)  78.7 W28-W30 N15-N20
      8 1667 1068 (N16W37 L= 86)   267     52        1 (1650/1038-1679/1100) -76.0 W35-W38 N13-N17
      9 1323 1049 (N04W16 L= 64)   208     89        0 (1303/1008-1356/1081)  80.5 W14-W16 N00-N06
     10  430  915 (S28E35 L= 13)   640    115        2 ( 377/ 886- 479/ 961)  40.6 E39-E31 S30-S24
     11  968  918 (S14E00 L= 48)   245    107        0 ( 922/ 908-1029/ 933)  21.8 E02-W03 S15-S12
     12  918  748 (S24W02 L= 50)   440    202        0 ( 824/ 732-1002/ 768)  21.8 E03-W07 S27-S21
     13 1215  703 (S18W21 L= 69)   859    186        1 (1132/ 657-1304/ 746)  52.1 W16-W24 S22-S13
     14  660  697 (S35E13 L= 35)   329    134        0 ( 598/ 688- 723/ 709)  26.6 E17-E08 S36-S32
     15  800  561 (S39W01 L= 49)   911    264        0 ( 710/ 492- 912/ 639)  61.9 E03-W06 S46-S30
     16 1129  582 (S27W21 L= 69)   213     62        0 (1107/ 552-1155/ 601)  76.0 W19-W21 S29-S25
     17 1151  235 (S48W48 L= 96)   572     69        0 (1122/ 200-1172/ 268) -14.0 W41-W53 S49-S45
     18 1297  181 (S44W68 L=117)  1058     87        1 (1281/ 135-1323/ 223)  -6.8 W59-W84 S46-S41
  

    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