Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  606 1574 (N24E41 L=167)  3295    574        6 ( 460/1447- 826/1748) -14.9 E70-E21 N19-N33
      2  977 1434 (N22E11 L=198)   327     99        1 ( 932/1415-1029/1462)   0.0 E14-E07 N21-N23
      3 1487 1349 (N28W24 L=233)   376    138        0 (1448/1296-1532/1398) -32.0 W20-W28 N25-N30
      4  832 1259 (N09E16 L=193)   796    300        1 ( 727/1214- 954/1339)  -3.6 E24-E08 N06-N11
      5  673 1186 (N02E24 L=185)   479     89        2 ( 626/1166- 727/1204)  18.4 E27-E20 N00-N03
      6 1129 1106 (N05W05 L=213)   484    208        0 (1104/1012-1151/1198) -79.7 W03-W06 N00-N11
      7  526  875 (S20E29 L=180)   641    136        0 ( 505/ 803- 553/ 933) -82.9 E30-E27 S23-S15
      8  602  755 (S26E22 L=187)   434    137        1 ( 538/ 727- 651/ 794)  -6.3 E27-E18 S26-S24
      9 1435  371 (S33W49 L=258)   701     99        2 (1381/ 354-1484/ 397)  12.5 W44-W53 S34-S30
     10 1290  365 (S37W37 L=246)   303     44        1 (1269/ 347-1315/ 389) -21.8 W35-W40 S37-S34
     11  822  189 (S68W10 L=219)   540     59        0 ( 800/ 164- 841/ 217) -31.0 W05-W16 S71-S64
  

    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