Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  770 1877 (N41E50 L=178)   925    251        1 ( 752/1829- 789/1929) -37.6 E57-E43 N37-N47
      2  598 1507 (N11E39 L=190)  1150    382        0 ( 539/1367- 680/1648) -84.3 E41-E36 N01-N21
      3 1525 1404 (N28W21 L=249)   271     36        0 (1496/1376-1548/1434) -21.8 W18-W23 N26-N29
      4 1376 1211 (N12W15 L=243)  1802    493        3 (1285/1004-1529/1431) -35.8 W04-W29 N04-N23
      5  755 1261 (S00E21 L=207)  1815    567        3 ( 538/1049- 863/1357) -18.4 E36-E09 S10-N05
      6  269 1272 (S09E57 L=171)   372     70        0 ( 246/1235- 298/1311)  68.2 E58-E54 S11-S06
      7 1824 1132 (N22W50 L=278)   823    118        2 (1798/1068-1840/1187) -50.2 W46-W51 N18-N25
      8  674 1064 (S13E21 L=207)   981    214        2 ( 576/1021- 779/1114)  39.3 E27-E15 S17-S08
      9 1819  890 (N08W56 L=284)  1340    206        0 (1793/ 815-1846/ 991) -52.1 W50-W59 N03-N13
     10 1053  697 (S25W11 L=239)  1015    392        0 ( 857/ 669-1246/ 731)  26.6 E01-W21 S30-S18
     11  505  697 (S39E27 L=201)  1615    422        5 ( 347/ 679- 679/ 724)  17.8 E41-E13 S43-S34
  

    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