Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  443 1467 (N19E48 L=211)   526    143        1 ( 417/1418- 464/1512) -39.8 E51-E44 N16-N21
      2  665 1447 (N21E30 L=228)   682    113        1 ( 642/1389- 694/1503) -81.9 E32-E28 N17-N24
      3  564 1303 (N10E34 L=225)   725    141        0 ( 528/1248- 608/1369) -41.2 E37-E29 N07-N14
      4 1341 1235 (N17W17 L=276)   271     98        1 (1288/1193-1370/1269) -26.6 W13-W19 N15-N19
      5 1063 1219 (N12E00 L=258)   543    182        1 ( 958/1192-1148/1251)   5.2 E07-W04 N11-N13
      6  273 1020 (S11E52 L=206)   789    173        1 ( 256/ 932- 290/1105) -78.7 E53-E50 S15-S05
      7 1166  968 (S01W09 L=268)  1191    298        1 (1097/ 841-1263/1102) -63.4 W04-W17 S07-N07
      8  627  922 (S12E23 L=235)   511    152        2 ( 568/ 904- 720/ 940)   5.7 E27-E17 S12-S10
      9  215  883 (S20E58 L=201)   474    147        0 ( 178/ 858- 255/ 926) -18.4 E62-E53 S21-S17
     10  903  868 (S11E05 L=254)   421    173        1 ( 865/ 800- 985/ 913)  45.0 E07-W00 S15-S07
     11 1239  702 (S16W18 L=277)   975    221        4 (1165/ 597-1308/ 788)  65.2 W15-W21 S23-S10
     12 1065  534 (S30W11 L=269)   525    118        2 (1007/ 511-1138/ 564)  18.4 W06-W15 S32-S27
  

    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