Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  858 1708 (N35E23 L=196)   638     59        0 ( 820/1684- 909/1732)  31.0 E25-E19 N33-N37
      2  959 1595 (N28E13 L=206)   291     69        0 ( 930/1562- 994/1617)  59.0 E13-E10 N25-N30
      3 1166 1399 (N18W04 L=223)   211     21        0 (1148/1378-1185/1419) -26.6 W02-W04 N16-N19
      4  761 1204 (N01E18 L=201)   656    253        0 ( 709/1094- 800/1328) -56.3 E22-E14 S04-N08
      5 1615 1212 (N14W35 L=254)   268     56        2 (1588/1189-1643/1239)   0.0 W32-W36 N12-N15
      6 1909 1143 (N16W63 L=282)  1612    181        3 (1862/1099-1947/1190) -33.7 W57-W68 N13-N19
      7  576  974 (S14E27 L=192)  1004    330        1 ( 465/ 865- 702/1089)  45.0 E34-E19 S21-S06
      8 1916 1055 (N11W65 L=284)   640    166        0 (1872/1036-1942/1084)  29.1 W59-W68 N08-N13
      9  370  985 (S15E42 L=177)   809    280        3 ( 267/ 933- 469/1033)  19.4 E51-E34 S18-S11
     10 1214  811 (S16W14 L=234)   455    169        2 (1149/ 755-1253/ 856) -29.7 W10-W17 S18-S13
     11  502  653 (S35E32 L=187)   401     61        1 ( 471/ 628- 530/ 680) -33.7 E34-E29 S36-S33
     12  655  452 (S47E21 L=198)   396    110        1 ( 600/ 433- 714/ 468)  10.3 E26-E15 S48-S45
     13  864  214 (S65W03 L=222)   503     80        1 ( 837/ 195- 892/ 247)  63.4 W00-W07 S67-S61
  

    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